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In order to investigate the nature of the antiferromagnetic structures in perovskite RMnOa, 
we study a Heisenberg J1-J2 model with bond alternation using analytical and numerical 
approaches. The magnetic phase diagram which includes incommensurate spiral states and 
commensurate coUinear states is reproduced. We discuss that the magnetic structure with spin 
ttii configuration (E-type structure) and the ferroelectricity emerge cooperatively to stabilize 
this phase. Magnetoelastic couplings are crucial to understand the magnetic and electric phase 
diagram of RMnOg. 

KEYWORDS: Heisenberg model, bond alternation, E-type structure, spin spiral state, A-type structure, 
JJMnOs, multiferroics 



Strongly correlated electron systems show various 
phase transitions involving complex order parameters, 
in general. Typically, in perovskite manganites -RMnOa, 
phases with spin, charge, orbital and lattice degrees of 
freedom emerges and competes with each other, through 
control of carrier concentrations and ionic radius of R 
ions.^ In these Mn ions, 3d electrons exhibit varieties of 
phases in nearly identical systems with a slight change 
in a few set of parameters. 

One of the recent topics in these /JMnOs is the pres- 
ence of multiferroic properties. For compounds with 
R =La, Nd, Sm, . . . , a spin A-type antiferromagnetic (A) 
phase is observed. In this phase, spins order ferromagnet- 
ically in the ab plane and antiferromagnetically along c 
axis, i.e., the magnetic ordering vector is q = (0, 0, 27r) in 
the orthorhombic lattice with GdFeOa-type distortions. 
However, when one substitutes R ions with smaller ones 
R = Tb, Dy, . . . , an incommensurate spin spiral (ICS) 
phase emerges accompanied by ferroelectricity.^'^ In this 
phase, magnetic propagation vectors are along the b axis, 
q = (0, Qm, 27r) where < q,„ < tt. The origin of the fer- 
roelectricity is the inverse Dzyaloshinsky-Moriya (DM) 
mechanism for the ICS structures. Here, a local elec- 
tric polarization is generated through the antisymmetric 
magnetoelectric (ME) coupling Pij cx e.^ x {Si x Sj), 
where Si is the spin direction at i-th site and e.^ denotes 
the unit vector connecting i and j sites. One of the ex- 
perimental evidences is that the reorientations of P from 
J* II a to i-* II c by chemical substitutions or application 
of magnetic fields are always accompanied by ab to be 
cycloidal-plane flops. '''^ 

Among various attempts to understand phase dia- 
grams of i?Mn03,''''^ finite temperature phase diagrams 
including ab- and 6c-cycloid phases have only been re- 
produced successfully by a model described as follows, 
so far: The model is a classical Heisenberg model defined 
on an three-dimensional orthorhombic lattice with near- 
est neighbor (n.n.) ferromagnetic exchanges Ji, and next 
nearest neighbor (n.n.n.) antiferromagnetic exchanges 



J2, as well as single ion anisotropics and DM interac- 
tions.^"' ^-'^ Models based on the J1-J2 Heisenberg model 
also explain spectra for magnons and electromagnons in 
these compounds. -'^^"^^ Construction of a realistic and ac- 
cessible microscopic model is quite appreciable from the 
viewpoint of comprehension of mechanisms as well as 
predictions and materials designs for novel phenomena. 

For smaller ions at i?=Ho, Tm, Lu, . . . , where J2/J1 
is expected to be larger, spins exhibit an E-type antifer- 
romagnetic (E) phase with spin ttii coUinear structure 
along the 6-axis with the magnetic propagation vector 
q — (0,7r, 27r). This phase also exhibits ferroelectricity 
P II a and thus is multiferroic. It is considered that 
the ferroelectricity is driven by the E phase orderings 
through the symmetric ME coupling \P\ cx {Si ■ Sj) of 
the magnetoelastic origin, ^^'^^ as illustrated in Fig. 1. 

The Heisenberg model with J1-J2 interactions alone, 
however, does not reproduce the E phase, as recently 
emphasized by Kaplan. A possible way to stabilize 
the E phase is to introduce a uniaxial anisotropy which 
enhances coUinear behaviors. The authors have studied 
the Heisenberg model with anisotropies and DM inter- 
actions, which successfully reproduces the A-ICS tran- 
sition, in the large J2/J1 region. Within the realistic 
range of parameters, however, we fail to observe the E 
phase. Another candidate is the biquadratic interaction 
which also enhances coUinearity.^^ At this point, how- 
ever, it is not clear weather such an interaction domi- 
nantly acts to stabilize the E phase in RMnOs. 

In this paper, we introduce an alternative model to 
investigate the spin structures of i?Mn03. Namely, we 
study a classical Heisenberg model with bond alterna- 
tions as a result of the magnetoelastic couplings men- 
tioned above. We clarify that the model can reproduce 
the A-ICS-E phase transitions in i?Mn03, and discuss 
the nature of the ME phase diagram. Since -RMnOa is a 
rare system which exhibits both symmetric and antisym- 
metric ME couplings E and ICS phases, respectively, it is 
also quite interesting to study the phase transition across 
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Fig. 1. Schematic view of the lattice distortion in E-type spin 
structure phase of RMnOs. Gray circle with arrow shows Mn 
ions with spin tt-U- structure. Displacement of oxygen ions 
(filled circle) gives rise to a uniform electric polarization. Due to 
the Goodenough-Kamamori mechanism, nearest neighbors ex- 
changes are modulated as ±AJi, which gives a gain in total 
energy for E-type phase. 




Fig. 2. Bond alternating J1-J2 model, depicted on (a) the ab 
plane and (b) the orthorhombic lattice. Thick and thin diago- 
nal lines represent nearest neighbor bonds with exchange cou- 
pling (1 -I- A)Ji and (1 — A)Ji, respectively. Dotted lines rep- 
resent bonds with J2. Interlayer antiferromagnetic coupling Jc 
is shown by dash-dotted lines, a, b and c represent orthorhom- 
bic axes. Thick arrows in (a) depict a spin configuration in the 
E-type structure. 



these phases in order to make further comprehension for 
the ME effects in these compounds. 

We study a classical Heisenberg model with spin ex- 
change bonds depicted in Fig. 2 as 

7^ = - ^ JijS, ■ Sj - «/2 ^ Si ■ Sj 
. ^ " 



<i,j> 



(1) 

<i,j> 

Here, the first summation is taken over n.n. bonds on 
the ab plane, and Jij takes the alternating values along 
6-axis in the form 



J^J = Jl(l±A), 



(2) 



as depicted in Fig. 2. Bond alternation parameter A is 
restricted within ^ A ^ 1. The second summation is 
for n.n.n. bonds along 6-axis. We take Ji > (ferromag- 
netic) and J2 < (antiferromagnetic). The third term 
of the Hamiltonian is the inter-plane antiferromagnetic 
couplings along c-axis. At T = 0, Jc merely create stag- 
gered stacking of a6-plane spin structure along c-axis, 
and thus irrelevant for the phase diagrams with respect 
to Ji and J2 . For example, it is determined automatically 
that a ferromagnetic alignment of spins in the ab plane 
stack antifcrromagnctically form the A-phase. Therefore, 
we may focus on spin structures within the ab plane. 

In the absence of the bond alternation A = 0, the 
model at r = gives an ordinary spiral state with On = 
nO, where 0„ describes the angle of the spin within the 
spiral plane. Here, n describes the position of the spin 



along the propagation vector || b. The rotation angle 9 is 
determined by cos 6* = 1/(27), where 



7 = -J2/J1 



(3) 



and the b component of the propagation vector is given 
by Qm = '2.9. 

At A 7^ 0, staggered modulation of Ji is introduced. 
Then wc introduce a variational spin state with uniform 
and staggered component for the rotation angle 



9n=ne+ [-1^5/2. 



(4) 



Total exchange energy per site scaled by Ji is described 
by e{9,5), which is calculated as 

e{9, 5) = -(1 + A) cos(e' - 5) - (1 - A) cos(6i + 5) 

7 cos 26* -7c, (5) 

where 7c = —Jc/Ji - Minimization of £{9, 5) with respect 
to 5 and 9 through de/dO = de/d6 ^ leads to 

(6) 



tan (5 ~ A tan 6', 



COS! 



'l__A2 

472 



A2 



1 



(7) 



From Eq. (6) we have (5 = in the limit 6 = 0, which 
implies an A-phase, irrespective of A. Similarly, at S = 
Tr/2 we have 9 = tt/2 which makes a collinear E phase. 
Critical value of A for the A-ICS boundary is given from 
Eq. (7) at 61 = as 

Al^) 



whereas at 9 



Since Ac^'' > 0, E phase is stabilized only in the presence 
of the bond alternation A. In Fig. 3 wc show the phase 
diagram. 

In order to justify the above analytical discussion, we 
perform the Monte-Carlo calculations at low tempera- 
tures. In the above, we have assumed the spin configu- 
ration of Eq. (4). On the other hand, the Monte-Carlo 
calculation does not restrict the spin configuration, and 
thus provides unbiased results. The calculation indeed 
confirms that only the three magnetic phases, i.e.. A, 
ICS, and E phases arc possible within our model so that 
the assumption of Eq. (4) is justified. As its consequence, 
the analytical results are precisely reproduced by the cal- 
culation. 

In our numerical calculations, we also add the "Hd = 
D^jS'^j, which makes magnetization along the z axis 
hard. Because of this term, the spins in the ICS phase 
rotate in the xy plane. This term makes the calculations 
stable by suppressing thermal fluctuations of the spiral 
plane without affecting the A-ICS-E transitions at T = 0. 
We take Jc/Ji=l and D/Ji^O.2 in the calculation. We 
analyze this model using the Monte-Carlo technique for 
systems with 48x48x6 sites with periodic boundaries. 

The calculation successfully reproduces the analyti- 
cally predicted phase diagram of Fig. 3. The circle, tri- 
angle and square symbols in Fig. 3 denote the points at 



v/r^, (8) 

7r/2 wc have the ICS-E boundary 



7- 



(9) 
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Fig. 3. (Color online) Phase diagram of the alternate J1-J2 model 
at the ground state. A, E, and ICS denote A-type, E-type, and 
incommensurate spiral phases, respectively. Solid lines are the 
phase boundaries determined by the analytical calculation. Cir- 
cle, triangle, and square symbols denote the points at which the 
A, E, and ICS spin structures are respectively obtained in the 
Monte-Carlo calculation at T / J\=0.2. 



which the A, E, and ICS spin structures are respectively 
obtained in the Monte-Carlo calculation at T/Ji=0.2. 

These magnetic structures are assigned from the 
peak position of the spin correlation functions at q = 
(0, 27r). We also measure the spin-hclicity correlation 
function to identity the phases. Here, the local spin- 
helicity vector is defined as hi = [Si x Si-^x + Si x 
Si+y)/2S'^ where x and y point the n.n. sites along [110] 
(pseudo-cubic x) and [-110] (pseudo-cubic y) directions, 
respectively. In spiral structures, the rotating spins give 
rise to a ferro-arrangement of the spin helicities, which 
results in a peak of the spin-heHcity correlation at q=0, 
while in coUinear phases, the spin-helicity correlation has 
no peak structure. 

Let us now show some details of the calculation. In 
Fig. 4, we show the temperature dependence of specific 
heat CsiT) = {l/N){d{'H))l{d{kB.T)), as well as the 
spin-helicity vector hs{T) where ha{T) {a—x, y and z) 
denotes the a component of the spin-helicity vector, at 
(a) (A, 7)=(0.4, 0.2), (b) (A, 7)=(0.4, 0.6), and (c) (A, 
7)=(0.4, 1.2) where the ground states are predicted to 

be (a) A-type, (b) ICS, and (c) E-type, respectively 

see also the phase diagram in Fig. 3. 

In all cases, we can see a single phase transition from 
paramagnetic (PM) phase at high temperatures to each 
ordered phase at low temperatures, at which Cs{T) ex- 
hibits a sharp peak. We also confirm that the choice of 
T/Ji = 0.2 gives a sufficiently low-temperature state. 

In the present calculation, the spins in the ICS phase 
rotate in the xy plane since we incorporate the hard-axis 
type spin anisotropy along the z axis in the Hamiltonian. 
Under this circumstance, h^iT) has a large value, while 
hx{T) and hy{T) are almost zero. We can indeed see that 
hz{T) starts increasing at the transition to the ICS phase 
in Fig. 4(b). On the other hand, as shown in Figs. 4(a) 
and 4(c), z component of the helicity is as small as other 
components in the coUinear spin phases. 

In Fig. 5, we show the calculated 7 dependence of the 
momentum qm/T^ as well as the spin-helicity for (a) 
A=0.6, (b) A=0.4, and (c) A=0.2. The data for the 
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Fig. 4. (Color online) Calculated temperature profiles of specific 
heat Cs(T) and spin-helicity vector hs{T) for typical points in 
the phase diagram where the ground states are estimated to be 
(a) A-type, (b) ICS, and (c) E-type. PM denotes the paramag- 
netic phase. 




Fig. 5. (Color online) Calculated 7 dependence of the momen- 
tum 5m /i" at which the spin correlation function has a peak at 
T/Ji=0.2, and 7 dependence of the spin-helicity z component 
at r=0 for (a) A=0.6, (b) A=0.4, and (c) A=0.2. The data 
of at r=0 are obtained by extrapolating the calculated tem- 
perature profile of hz{T). Solid and dashed lines are qm/i" and 
hz obtained from the analytical approach. 
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spin-hclicity arc obtained by extrapolating the calculated 
temperature profile of hz{T) to T/Ji — )• 0. Those data as 
functions of 7 consistently exhibit A-ICS-E phase tran- 
sitions, at least qualitatively. 

Furthermore, let us compare these data with the varia- 
tional spin structure given in Eq. (4) . The spin structure 
gives qm = 29 and = sin cos (5, where 9 and 5 are 
determined by Eqs. (6) and (7). As plotted in Fig. 5, 
the Monte-Carlo data for and agree well with the 
analytical results given by solid and dashed lines, respec- 
tively. 

So far wc have confirmed that the spin structure given 
in Eq. (4) as well as the A-ICS-E phase transition derived 
from that is quite valid. Then we proceed to consider 
the case that the bond alternation is driven by lattice 
distortions to minimize the total energy. We assume that 
a lattice distortion r creates the bond alternation A = 
with g being the coupling constant. Total energy is given 
by a sum of the spin exchange energy Eg (A, 7) and an 
elastic energy oc r^, or equivalently 



Etot 



A2 

E.(A,7) + ^. 



(10) 



Here, A is the dimensionless spin-lattice coupling con- 
stant. Hereafter, we refer to the model as Peierls J1-J2 
model. 

The spin exchange energy is derived by applying 
Eqs. (6) and (7) to Eq. (5). In the A phase at A < Ai^\ 
the minimized energy is derived from 9 = 5 = Q 'as 



e,(A,7) = -2 + 7-7, 
.(E) 



(11) 



while in the E phase at A > Ac'^'', we have 9 ^ 5 = ti /2 
so that 



e,(A,7) = -2A-7-7,. 
In the ICS phase at A^^^ < A < A^^^ , we have 



es(A,7) = - 



1- A^ 
27 



+ 7 



1 + A^ 
1- A2 



7c- 



(12) 



(13) 



Bond alternation A is determined by minimizing etot 
with respect to A within the range ^ A ^ 1. The result 
in the region 7 > 1/2 can be summarized as follows. In 
Fig. 6, we schematically depict etot as a function of A for 
various A. We observe a transition from ICS phase to E 
phase at a critical coupling Ac given by 

Ac' = 1/(27). (14) 

Here, etot have degenerate energy minimums at A = 
(ICS region) and at A > A, (E region). As A is changed 
across Ac, we have a jump in the value of A for the en- 
ergy minimum and thus the transition is first ordered. 
Furthermore, there exists a characteristic coupling 



Ai' = vV+T -7, 



(15) 



where E phase has a local minimum in its energy with 
respect to A at A > Ai. This implies an existence of a 
metastable E state at A > Ai. Similarly, we also have 



As" = 27/(472 



1), (16) 
where an ICS state at A = is metastable at A < A2. 




Fig. 6. Phase diagram of the Peierls J1-J2 model. The solid curve 
shows the first-order phase boundary between E phase and ICS 
(F) phase at 7 > 1/2 (7 < 1/2). The vertical solid line at 7 = 
1/2 shows the second-order phase boundary between ICS and F 
phases. Hatched area shows the coexistence region. Inset: Total 
energy etot as a function of A plotted in arbitrary units, for the 
case of 7 > 1/2. Dimensionless coupling constant A is increased 
from top to bottom, as A < Ai, A = Ai, A = Ac, A = A2 and 
A > A2, respectively. 



Therefore, in the case 7 > 1/2, we may have coexistence 
of the ICS and the E phases at Ai < A < A2. 

Similarly, in the case 7 ^ 1/2, 2nd order transition 
between the A and the ICS phases is observed at Ac, 
while the metastable E state exist at A > Ai. The A 
states is always metastable at A > Ac in this case. The 
forms for Ac and Ai coincide with those for the previous 
case given in Eqs. (14) and (15), respectively. The phase 
diagram is summarized in Fig. 6. 

Let us now compare these results with the magneto- 
electric properties in i^MnOa. In these compounds, com- 
mensurate E phase is always accompanied with ferroelec- 
tricity, which is consistent with our result that E phase is 
stabilized only at A 7^ 0. Although it has previously been 
discussed that ferroelectricity is triggered by magnetic 
E phase through lattice distortions, our result suggests 
that the magnetic E phase and ferroelectricity coopera- 
tively emerge to stabilize themselves. Therefore, Peierls- 
type spin-lattice couplings are essentially important to 
understand the electric as well as the magnetic phase di- 
agram of i^MnOa. Recent report shows that, in i?Mn03 
there exists a region of possible coexistence for ICS and E 
phases. The Peierls model discussed here indeed shows 
such a coexistence at around the Ist-order transition 
points between ICS and E phases. Experimental data 
in further details should give crucial tests for our present 
results. Theoretical approaches based on more realistic 
models including anisotropics, DM and biquadratic in- 
teractions should also be important to understand the 
whole phase diagrams of i^MnOa in detail. 

The author would like to thank S. Miyahara for stim- 
ulating discussions, as well as S. Ishiwata, Y. Taguchi 
and Y. Tokura for suggestions from experimental points 
of view. This work was partially supported by Grant-in- 
Aid for Scientific Research as well as High Tech Research 
Center Project from the MEXT, Japan. 



J. Phys. Soc. Jpn. 



Letter 



N. FURUKAWA, M. MOCHIZUKI 5 



1) For a review, "Colossal Magnetoresistive Oxides", ed. by Y. 
Tokura (Gordon & Breach Science Publisher, 2000). 

2) T. Kimura, T. Goto, H. Shintani, K. Ishizaka, T. Arima and 
Y. Tokura: Nature 426 (2003) 55. 

3) T. Kimura, S. Ishihara, H. Shintani, T. Arima, K. T. Takahashi, 
K. Ishizaka, and Y. Tokura: Phys. Rev. B 68 (2003) 060403. 

4) H. Katsura, N. Nagaosa, and A. V. Balatsky: Phys. Rev. Lett. 
95 (2005) 057205. 

5) L A. Sergienko and E. Dagotto: Phys. Rev. B 73 (2006) 094434. 

6) M. Mostovoy: Phys. Rev. Lett. 96 (2006) 067601. 

7) Y. Yamasaki, H. Sagayama, T. Goto, M. Matsuura, K. Hirota, 
T. Arima, and Y. Tokura: Phys. Rev. Lett. 98 (2007) 147204. 

8) Y. Yamasaki, H. Sagayama, N. Abe, T. Arima, K. Sasai, M. 
Matsuura, K. Hirota, D. Okuyama, Y. Noda, and Y. Tokura: 
Phys. Rev. Lett. 101 (2008) 097204. 

9) S. Dong, R. Yu, S. Yunoki, J.-M. Liu, and E. Dagotto: Phys. 
Rev. B 78 (2008) 155121. 



10) M. Mochizuki and N. Furukawa: J. Phys. Soc. Jpn. 78 (2009) 
053704. 

11) M. Mochizuki and N. Furukawa: Phys. Rev. B 80 (2009) 
134416. 

12) S. Miyahara and N. Furukawa: preprint, arXiv:0811.4082. 

13) R. Valdes Aguilar, M. Mostovoy, A. B. Sushkov, C. L. Zhang, 
Y. J. Choi, S-W. Chcong, and H. D. Drew: Phys. Rev. Lett. 
102 (2009) 047203. 

14) M. Mochizuki, N. Furukawa and N. Nagaosa: arXiv:1001.3905. 

15) T. Arima, A. Tokunaga, T. Goto, H. Kimura, Y. Noda, and 
Y. Tokura: Phys. Rev. Lett. 96 (2006) 097202. 

16) Ivan A. Sergienko, Cengiz Sen and Elbio Dagotto: Phys. Rev. 
Lett. 97 (2006) 227204. 

17) T.A. Kaplan: Phys. Rev. B80 (2009) 012407. 

18) S. Ishiwata, Y. Kaneko, Y. Tokunaga, Y. Taguchi, T. Arima 
and Y. Tokura: arXiv:0911.4190. 



